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I. INTRODUCTION 

Studying Quantum chromodynamics (QCD) phase transition and properties of hot/dense 
quark matter at high temperature and baryon density is one of the most important topics 
of high energy nuclear physics. The Relativistic Heavy Ion collider (RHIC) and the Large 
Hadron Collider (LHC) provide the opportunity to investigate properties of nuclear matter 
at high temperature and small baryon density. It is now believed that the system created 
at RHIC/LHC is a strongly coupled quark-gluon plasma (sQGP) and behaves like a nearly 
“perfect” fluid [1, 2], One crucial quantity is the shear viscosity over entropy density 77 /s, 
which is required to be very small to fit the elliptic flow at RHIC/LHC. The result from 
AdS/CFT correspondence gives the lower bound of 77 /s = -A [3, 4], which is very close to 
the value used to fit the elliptic flow v 2 [5-7]. 

The anti-de Sitter/conformal field theory (AdS/CFT) duality [8-10] is discovered through 
pioneering works on the near horizon structure of black branes (or black holes) [11-17] 
and the scattering process of branes and bulk probe fields [18-21], It is generalized to 
nonconformal brane background in the near horizon limit in ref. [22], which is called the 
gauge/gravity duality nowadays. The gravity/gauge duality or AdS/CFT correspondence 
provides a revolutionary method to tackle the problem of strongly coupled gauge theories. 
It has been widely used to investigate QCD phenomenology, e.g. glueballs [23-25], hadron 
spectra [26-29, 31], deconfinement phase transition [27-30] and transport properties [35]. 

The shear viscosity in AdS/CFT was firstly calculated in Ref. [35] through relations 
between Green-Kubo formula 1 and the absorption cross section of gravitons [18-21]. Studies 
on the near equilibrium QGP from AdS/CFT duality in [36] gives a recipe of extracting two- 
points real-time thermal correlators via classical bulk action. Following [36], the authors of 
[39, 40] calculated the first order transport coefficients in near extremal D3 brane background 
and found that, in the long-distance and low-frequence limit, these correlators turn into 
hydro dynamical forms. Second order transport coefficients of this system were calculated in 
[41]. The framework that investigates transport properties of a fluid via its corresponding 
gravity is called the fluid/gravity correspondence, and the most notable feature in above 
works [39-41] is the use of Green-Kubo formula, thus one may call it the Green-Kubo 
formalism of fluid/gravity correspondence. 

While the Green-Kubo formalism becomes popular in extracting transport properties of 
liquid-like plasma 2 , another systematical and powerful formalism—the boundary derivative 

1 which is first proposed by Ryogo Kubo [37] in statistical mechanics and recast into field theory formalism 
by Akio Hosaya et.al. in [38]. 

2 the literature on this topic can be found in the references of [42]. 
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expansion (BDE) formalism [43, 44] has been developed. The most remarkable feature for 
this formalism is the use of boundary dependent boost parameters for the bulk metric in 
the (in-going) Eddington-Finkelstein coordinate. Expansions are implemented with respect 
to boundary derivatives of the boost parameters and all the dissipative terms of boundary 
fluid are metric perturbations (in the large r limit) solved from Einstein equation. The first 
example of the BDE formalism of fluid/gravity correspondence is a duality of AdS§ black 
hole in the bulk to finite temperature conformal J\f = 4 SYM plasma in the boundary [43], 
where the transport coefficients were calculated to second order. The BDE formalism was 
applied in several other models with AdS§ black hole background: 1. The kldS^-dilaton 
model [45], where the gravity side is an AdS§ black hole plus a boundary dependent dilaton 
field, while in the boundary is a fluid with forcing terms. 2 . The charged AdS 5 black hole 
model [46, 47], where the bulk is a charged AdS$ black hole, and the fluid on the boundary 
has a chemical potential. The addition of the Chern-Simons term for the £/(1) gauge field 
cause the appearance of vorticity in the first order dissipative expansion of R-charge current. 

The development of fluid/gravity correspondence interweaves with the studies on hydro- 
dynamical modes on the world-volume of the blackfold [48-51], which opens a window to 
extract the dynamical information on the world-volume of black branes in flat spacetime. 
The most obvious difference of this kind of research from the fluid/gravity correspondence 
is the need of Dirichlet boundary condition on a hnite cut-off surface. In [52], the effective 
hydrodynamics on ap + n + 2 dimensional “rigid wall” located atr — RinaD — p + n + 3 
dimensional spacetime is studied, where p and n + 1 are the number of (spatial) dimension 
of the brane and the sphere, respectively. Based on this, [53] studies the effective hydrody¬ 
namics on the world-volume of black D3 brane to first order, both the thermodynamical and 
the viscous quantities depend on the location of the cavity and the horizon, however r\js of 
this model is still -A. Ref. [54] investigates the effective hydrodynamics of rotating black 
D3 branes. The common feature of [53, 54] is the use of dimensional reduction, which trans¬ 
forms the effects of transverse directions into massless holds on the longitudinal directions 
in which the branes lie. Thus makes us focus on the worldvolume theory. 

The Green-Kubo formalism and BDE formalism of fluid/gravity correspondence provide 
powerful tools for us to study the liquid-like QGP (see e.g., [55, 56] for a phenomenological 
review on heavy ion collisions for theorists). Generally speaking, QGP is a liquid-like plasma 
with small shear viscosity, and can be described by relativistic hydrodynamics quite well. 
Results from lattice show that QGP exhibits nonconformal properties especially around the 
critical temperature T c , e.g. the shear viscosity over entropy density ratio has a minimum 
around T c , and the bulk viscosity over entropy density shows a peak around T c [57-60]. 
This behavior has been described in bottom-up holographic QCD models [61-64], However, 
current studies using fluid/gravity duality from top-down method are mostly on AdS back¬ 
ground, whose dual fluid is of course conformal, thus may only reflect the properties of QGP 
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at the conformal regime, i.e., above 2 T c . However, when we are at the nonconformal zone 
around T c , AdS gravity may no longer be proper for a nonconformal gauge theory. 

A natural choice in top-down holographic way to tackle the strongly coupled nonconformal 
plasma is to build models using the nonconformal D-branes. Such studies include e.g. [65] 
for Dl-brane and [66] for Dp-branes with p > 2, where the Green-Kubo method is used. 
There is another interesting work on this topic which can handle more cases including 
p = 0,1 and the fundamental strings in type II string theory (but p = 5 excluded): the Ref. 
[67], where the BDE formalism in Fefferman-Graham coordinates developed in [68] is used. 
Besides the Dp-brane backgrounds like in Refs. [66, 67], one may also use the compactihed 
D-brane backgrounds, e.g. the compactihed D4-brane. The compactihed D4-brane is the 
background of D4-D6 model [29] and the Sakai-Sugimoto (SS) model [31], which are two 
nonconformal holographic QCD models from top-down. The background of these two models 
is the compactihed D4 black branes under the near horizon limit. There are some previous 
studies on the transport properties of this background. The sound speed and bulk to shear 
ratio was calculated in [69], the shear to entropy density ratio was argued in this reference 
to be 1/47T by showing that the SS model background is in the class of [71]. Using the null 
horizon focusing equation, Ref. [70] also calculated the ratio of the bulk viscosity to the 
shear viscosity. 

Based on the above review of the relevant literatures, one can see that there is a lack 
of parallel formulation with [43]. This motivates us to develop the framework of BDE 
formalism of fluid/gravity duality for nonconformal gauge theory plasma. In this paper, 
by using the BDE formalism of fluid/gravity correspondence, we calculate the first order 
transport coefficients of the nonconformal QGP under the quenched limit on the boundary 
of both the D4-D6 and SS model’s background. Our results are consistent with former 
studies by other methods. The previous results together with ours reveal that the plasma 
of the D4-D6 and SS model is nonconformal with a small bulk viscosity and saturates the 
KSS bound [3, 4], and this agreement guarantees the validity of the BDE formalism of 
fluid/gravity correspondence for nonconformal D-brane backgrounds with more than one 
submanifolds reduced. This work can been seen as a nonconformal counterpart that is 
parallel with Bhattacharyya et al.’s AdS§ construction in fluid/gravity correspondence [43]. 

This paper is organized as follows: After introduction, we will firstly give the preliminaries 
from 10D compactihed black D4 brane background to a 5 dimensional one in section 2 in 
order to make connection with the recipe of fluid/gravity correspondence. Then, in section 
3, we will solve all the first order perturbative ansatz and get the metric which perturbatively 
solves the Einstein equation to the first order. By making use of this solution, we calculate 
the boundary stress tensor for the QGP which corresponds to the bulk of SS model in section 
4 and analyze its transport properties. We give the discussion and outlook in section 5. 
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II. THE SETUP 


In this section, following [69], we will show how to derive the action and classical back¬ 
ground of the D4-D6 and SS model into 5D form through dimensional reduction as in 
[53, 54], The purpose of doing this is to make connection with [43], more details can be 
found in Appendix A. 

The D4 brane action of type IIA theory in Einstein frame is given as 


5 = 


?K 2 

0 


d w xV^G 


R (io) _ I(io V 0)2 _ ~ 2 

2 ’ 2-4! 


iov7^2 9s e l p 2 


( 1 ) 


where 2 k^ 0 = (27i) 7 g 2 l 3 is the 10D gravitational coupling and 10 V stands for 10D nabla. G 
is the determinant of the following diagonal 10D metric tensor: 


ds 2 = e 2aiA g MN dx M dx N + e 2aaA (e 2PlB dy 2 + e 2 ^ B ^ ab dO a dO b ) , (2) 

where gMN, A and B only depend on x M , the coordinates of first 5 dimensions, and y a b with 
a, b = 1, 2 , 3,4 is the metric on the S 4 . and j3\^ are four parameters whose value will be 
clear in the following context. The explicit form of dilaton and RR field are given in (19). 
It should be noticed here that y is also a compact dimension and we will integrate out both 
y and the 4-sphere to get a 5D effective theory. 

From Eq.(2), we have y/—G = e 5 (“i+« 2 )A+(/ 3 i+ 4 fc)B w pj 1 ^ _ dety^ the deter¬ 
minant of the metric on unit 4 sphere. During the reduction process, we have used the 
following relation: 

S ~ J d 10 xV^G(R iw) + •••) = I d 5 xy/^e {3ai+5a2)A+ ^+ 4 M B (R + •••). (3) 

To avoid the appearance of non-minimal coupling of the gravity with the scalar field in the 
reduced theory, one should set 

5 

an = — g, «2 = 1, fd i = 4, /?2 = —1, (4) 

so Eq.(2) becomes 

ds 2 = e~^ A g M Ndx M dx N + e 2A+8B dy 2 + e 2A ~ 2B dVL\. (5) 


From Eq.(A7), the 10D Ricci scalar has the form of: 


R m = e f A 


R + — V 2 A - — (dA) 2 - 20 (dB) 2 
3 3 


+ 12e 


-2A+2B 


( 6 ) 
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During the reduction process, we have 




( 7 ) 


= \/—< 7^/7 (r+ y V 2 yl - y(cb4 ) 2 - 20 (8B) 2 + 12e~f A+2 A , ( 8 ) 


\Z=G( 1 O V0 ) 2 = y/^-g^e-f A G MN d M( t>d N ^ = /^^W) 2 , 




Ql ^-^A+8B 


(9) 

( 10 ) 


Therefore the D4 brane action Eq.(l) is reduced to its 5D form and takes the form as 


5 = 


2«§ 


J d 5 x^g 


R - i(«90 ) 2 - y {dAf - 20 (dB) 2 - V(<j>, A, B) 


V(<f>,A,B) = %1-fA+sB _ 12e -fA+2B, 


( 11 ) 


where fi ) 5 is the 5D surface gravity with the following definition: 


1 _ D1D4 

24 = 24 / 


( 12 ) 


with V\ — f dy the volume of the compact circle. The system turns into a 5D Einstein 
gravity minimally coupled with 3 scalars (j>,A,B, and V(<p,A,B) is the scalar potential. 
The EOMs for this reduced system are: 


E 


MN 


T M N — 0, 


v 2 0 - A+8B = 0 , 

4 


V 2 kL + 


17 Ql C i-^A+8B _ ]^ e -fA+2B _ Q 
80 5 


Vie - 9l e %-fA+8B + 3 e -fA+2B = 0j 
10 5 


( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 


where 


Emn = Rmn — ^ m nR 


( 17 ) 


is the Einstein tensor in the 5D spacetime, and 


Tmn = - ydM<pd]\r(f> — -5 , mv(^ 0) 2 ^ + y ^SmABnA — -gMN(dA) 2 
20 ( DuBOnB — -gMN(dB) 2 j — -gMNV, 


( 18 ) 


which can be viewed as the energy-momentum tensor in the bulk. 
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The classical solution for black D4 brane in Einstein frame reads 

3 5 Ar 2 _3 5 

ds 2 = H 4 *(-f(r)dt 2 + dx 2 ) + Hf— + H 4 *dy 2 + Hjr 2 d^l 

fir) 

3 ^ 

e * = e *-*o = F 4 = ^-ig 4 e 4j H, = 1 + /(r) = 1 - If, (19) 

where = e $0 and Q 4 = (2nl s ) 3 g s N c /fl4. 3 Note that we write one of the directions that 
the D4 brane lies (denoted by y ) separatedly from the other 3 directions (denoted by {if}) 
in order to compare with (5). Under the near horizon limit, the above metric becomes 

dit z - (jj h i-firyU ' 1 + dx z j + (1) ^—J t ^L^ a dy 2 + L 1 ^Acl£l\, (20) 

e Mi) ! - < 21 > 

where L s = Q^/3 = ng s N c R s . The above metric differs from the D4-D6 model and the SS 
model for the interchange of t with y. Besides, it is in the Einstein frame, not string frame. 
Comparing Eq.(5) with Eq.(20), we have 

e^ = Lsorso 5 e =L “rw, (22) 


and the reduced 5D metric is 

T 4 

ds 2 = Lri (—/ (r)dt 2 + dx 2 ) H — 4 - dr 2 . (23) 

rs/(r) 


From its Ricci scalar R ~ ~ 6 r n / 3 (14r 3 + r^), when r — > 0, R will become minus infinity so 
r = 0 is the curvature singularity and away from that point the above metric will always 
be regular, thus we will only focus on the regime of r > 0 from now on. At the boundary 
r —> oo, R —> 0 so Eq.(23) is asymptotically flat, which is not obvious for the appearance of 
r 5 / 3 j n first 4 dimensions. 

We turn to the in-going Eddington-Finkelstein coordinate by making the following trans- 

r 3/2 

formation as dt = dv — then the above metric becomes 

ds 2 = Lr^ (—f (r)dv 2 + dx 2 ) + 2 L^r^dvdr. (24) 


r = 0 is still the curvature singularity of this 5D metric but everywhere away from that is 
regular. Since we have already lost the track of dimensions in the process of dimensional 
reduction Eq.(2), keeping L explicit will be insignificant, so from now on we set L — 1, 


3 The normalization condition for Q4 here is 2n 2 y4N c = f s4 F4, where 2 k 2 = 2 Kf 0 g s 2 and y 4 = ((2-7r) 4 Z 3 ) 1 
is the D 4 brane charge. 








which means Q 4 = 3. After a boost of coordinate: dv = —u^dx^, dx l = P'^dx 11 , where 
Ppv = Vi>M + we have 


ds 2 = r 3 (—f (r)u ll u u dx p dx u + P^dx^dx”) — 2r^u^ l dx tl dr, 


^ = 7(!,&), 7 


1 


(25) 


In the above metric, is the four-speed of the relativistic fluid with the normalization 
u = —1. P^ u is the projection tensor of the boundary with P /w P up = P p , which projects 
any tensor to the plane orthogonal to vJ l . As one can check, Eq.(25) is the Oth order 
solution of 5D EOM. The boundary of Eq.(25) is actually a fluid with constant temperature 
and velocity, which is of course in global equilibrium. 

In order to mimic slightly deviations from local equilibration and the anisotropy of the 
fluid, we promote the four parameters in Eq.(25) to be x p dependent: r# —> rn(x), u^ —> 
Ufj,(x), with the requirement that |^| -C 1, where T is the local temperature of the fluid. 
Then Eq.(25) becomes 


ds 2 = r 3 (— f(r H (x), r)u fl (x)u„(x)dx p dx l/ + P fiu (x)dx p ‘dx u ) — 2 reu fi (x)dx tl dr, (26) 


which is no longer the solution of 5D EOM, but we can make it as the solution again by 
putting some perturbations in. Using the method of [43], we should firstly expand the fluid 
quantities of Eq.(26) at some special point, say, x p = 0 in the local rest frame of the fluid 
as: 


tt/i = + x v d v fi j V ll , r H (x fl ) = r H ( 0) + x p d^r H . (27) 

r#(x M = 0) is the location of event horizon corresponding to x p = 0 in the boundary, it 
relates with the local equilibrium temperature of the fluid at that point. In order to keep 
the formulations neatly, we will just denote it as r# in the following calculations but one 
should always remember that it is a local quantity at x p = 0. Then we have 

u^dx p = —dv + x^d^ftidx 1 , u fJ ,u u dx p, dx L/ = dv 2 — 2x p d ll (3idx l dv, 

3r 2 

P lw dx ,l dx 1 ' = dx 2 - 2x p d ll /3 i dx l dv, f(r H (x),r) = f(r )- ^-x p d^r H . (28) 

Thus Eq.(26) becomes 

( 5 \ 1 

—r 3 / -1- T-x^d^rn I dv 2 - ^-x^d^idvdx 1 + 2 r&dvdr 

r 3 ) r 3 

+ Adx 2 — 2r^x p, d fJ ,(3 i dx l dr. (29) 

The above metric deviates the solution of Einstein equation slightly by the first order bound¬ 
ary derivatives at x p , we will see that adding some perturbation terms will make it the 
solution again, and these perturbations are solved in the next section. 
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III. THE FIRST ORDER PERTURBATIONS 


The SO(3) symmetry in Eq.(23) separates the perturbations into tensor, vector and 
scalars of SO( 3), and we will make use of this advantage to solve these three kinds of 
perturbations one by one. Generally speaking, all the perturbation ansatz will have the 
form: 


( 


difii, 


P(r) x d v p h 


for the scalar part; 
for the vector part; 


<7jj, for the tensor part, 


(30) 


where P{r) is some function of r and can be solved through Einstein equation with the 
boundary conditions as 


• P(r) is regular at r = r#; 

• lirn^oo ^ -A 0. 

Here n = 0 or n = 3 depends on the nature of perturbation terms. We can see that the 
perturbations will always be of the form of Eq.(30) with the above boundary condition 
implemented. 


A. The tensor part 


We set the tensor part perturbation as 

d s li)T = r^aij(r)dx l dx j . 


(31) 


The EOM that a tJ satisfies is: 

Eij - = T tj - ^8ij5 kl T kl . 

It turns out that the differential equation for cpj is 


d_ 

dr 



-5r2aij, 


(32) 


(33) 


where a %3 = <%/%) — \8ijd k Pk is the spatial part of shear stress tensor. The purpose for 
writing the EOM for tensor part like this is due to the traceless of the trace part of EOM 
should be removed from the diagonal components. The equation for the 1st order tensor 
perturbations takes similar form for different models, as can be seen e.g. from [43, 53, 54], 
the reason for this may be due to the universality of the shear viscosity in supergravity 4 


4 Tensor part perturbation corresponds to the shear viscous term. 
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[3, 71]. Since the metric of SS model is also in the class of [71], so it is natural for Eq. (33) 
to take such a form. We write a i3 as a t j = F(r)aij, then F(r) can be solved from 


F" 




H 77 / 


5 


r 4f 


F' + ~7 = 0 , 

r 2 J 


(34) 


from which the result can be solved as 

1 


F(r) = C 2 + 


3y/rH 


2V3 ( arctan 1 - hME _ Mctan 1 + 2 Eh~H 




V3 


In 


{y/r + y/rH) 2 (r + yfrrjj + r H ) 
(r - y/rr]l + r H ) 


+ C, In ■ 


H 


73 

- \n(y/r - 


(35) 


Regularity at r = t# requires C\ — 2 and the normalizability at r —> oo requires C 2 = 
2ir/^3rH 1 thus 


F(r) = 


3 y/r^ 


273 ^arctan 


1 - 2 y/r/r H 1 + 2 a Jr/r H 

— arctan- — -h n 


73 


73 


In 


(y/r + v^X) 4 ( r + y/ rr H + r H ) 2 (r 2 + rr H + r 2 H ) 


(36) 


It is regular at the whole regime of r > 0 and vanishes to 0 asymptotically. 


B. The vector part 


For the vector part, we set the perturbation ansatz as 

2 i 

osm ) v =- —Widx av. 

r 3 

The constraint equation for the vector perturbation is 


which gives 


g rv (E vi - T vi ) + g rr (E ri - T ri ) = 0, 


The dynamical equation is 


It turns out that wAr) satisfies 


diTH + 2 r H d v Pi = 0 . 


Eri ~ T ri = 0. 


2 5 r2 

w'! - w[- ^rd v fii 


l __ % 


r " 2r7 


= 0, 


H 


(37) 


(38) 


(39) 


(40) 


(41) 
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from which the solution is given as 

2 r ^ t 

w i( r ) =- 3 ~d v Pi + ^r 3 C u + C 2i . (42) 

r H ^ 

It is easy to see that the above general solution is regular at r H . The other boundary 
condition for the vector part perturbation is 

lim ^ 0, (43) 

which means Cu must be 0. The appearance of C 2 i will cause the (Oz) components of 
boundary stress tensor go out of Landau frame. So if one likes to express the boundary 
stress tensor in Landau frame, C- 2l should be set to 0. Thus the final result for the vector 
perturbation of 1 st order is 

dsfi)v = 4r® d v /3idvdx l . (44) 


C. The scalar part 


The scalar part, similar to other works on the effective hydrodynamics of black branes, 
e.g. [53, 54], is the most complicated part. We set the scalar part perturbation as 


dS(i) S 


k(r ) 

4 
7 * 3 


dv 2 + r3h(r)Sijdx l dx J + 2r&j{r)dvdr. 


(45) 


In our case, the gauge condition tr^^jr^p] = 0 [43, 46] cannot be used here for solving 
the scalar part perturbation, since this will cause inconsistence when solving the EOMs 
and make the surface stress tensor unrenormalizable. Other gauge condition like h(r) = 1 
[47] cannot be used either, since the spatial trace part of the metric is non-trivial in the 
nonconformal case here. Thus we need to keep all the 3 unknowns. However the labor cost 
by solving all of them gives us a bonus that there will be a bulk viscous term appeared in 
the surface stress tensor, which doesn’t appear in the conformal models with AdS gravity 
like in Refs. [43, 46, 47]. We have two constrain equations for the scalar sector: 


g rr (E rv - T rv ) + g rv (E vv - T vv ) = 0, (46) 

g rr (E rr - T rr ) + g rv (E rv - T rv ) = 0, (47) 

which separately give 

d v r H = —\ r H di(3i (48) 

5 

and 

3(5r 3 - 2 r%)ti - 30 r 2 j - 5 k' + = 0. (49) 
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We also have a total number of 7 dynamical equations for scalar perturbations, four of them 


3 

£ 

i —1 


^ 7 * 7 * 'J~ FT 0 ^ 

Erv T r v 0 } 

Evv Tyv O7 

{Eu - Tu) — 0 


(50) 

(51) 

(52) 

(53) 


come from the Einstein equation Eq-(13) and three of them 

V 2 0 - 9 e f-f -A+SB = g 
4 


V 2 A + 


^3- e ±-34 A+8B _ V2^_™ A+ 2B _ g 

5 


80 

V 2 £ - 


^-^o4+8S ^ _ifi 


+ -e-^ A+2B = 0 
5 


(54) 

(55) 

(56) 


are from the 3 scalar field equations in the 5D bulk, namely, Eqs.(14,15,16). This looks 
horrible at a first glance, but fortunately, not all of them give useful message. It turns out 
that Eq.(51) and Eq.(52) come out of linear compositions of specific constraints with Einstein 
equation of the scalar sector, so they are not independent equations, and Eqs.(54,55,56) give 
the same differential equation for the 3 unknown scalar perturbations. So we only need to 
solve Eqs.(49,50,53) and Eq.(54) to nail Eq.(45) down, among which the last 3 equations 


are 


0 = Qrh" + 9h' — 10/, 

0 = 12 r 4 fh" + 12(4r 3 - r 3 H )ti - Qrk" - 3k' 
— 6(5 r 3 — 2 r%)j' — 90r 2 j + 20 r^di/3 i: 

0 = 2 r 3 fj' + 12r 2 j + 2k' — 3 r 3 fh' — 2 


(57) 

(58) 

(59) 


We will choose Eqs.(49,57,59) to solve the 3 unknown scalar perturbations that we set in 
Eq.(45). From Eq.(49) we have 


6 r 2 j + k' = -(5r 3 — 2 r 3 H )h' + 2r^diPi, 
5 


(60) 


and after putting it into Eq.(59), we get 

—10r 3 /j' = (15r 3 + 3 r 3 H )h' + 10 r^d^i. (61) 

Then, putting the above equation into Eq.(57) one can finally get the equation for h\ 

d ( d „dh\ 5 3 _ „ 

T r f~T + o r2 ^A = 0 . 
dr \ dr J 3 


( 62 ) 
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Without losing of generality, we set h = F h (r)dif3i, and F h (r) satisfies 


d_ 

dr 


r 4 f 


dFh 

dr 


3 r2 ‘ 


(63) 


If comparing the above equation with Eq.(33), one can get F h = F/3 without solving it, 
thus 


h = -F(r)dA- 

Inserting h into Eq.(57) one has the equation for j = Fj(r)di/3i as 

10F' = 2 rF" + 3 F'. 


(64) 


(65) 


This is a 1st order differential equation, and the solution can be obtained by direct integra¬ 
tion, the result is 

1 T 1 — 2\/r Iru 

arctan ■ 


Fi(r) = - - - r -t + Q + 


5 r 3 - rjj 


30 yjr^ 


273 f arctan - _ 2 _ 1 + 2>/^ 




73 


73 


In 


(77 + - v /fy7) 4 (r + y/rr H + r H ) 2 (r 2 + rr H + r 2 H ) 


( 66 ) 


Since the above expression is already regular at r = r h, the remaining boundary condition 
for j is 


lim F; —y 0. 


(67) 


Thus Cj = 1 ^ 3 2_ , so we have finally 


Fj(r) = - 


2 rz — rh 1 ^ 

-— H- F 

5 H> -r% 10 


We set k = F k (r)di/3i likewise and substitute h and j into Eq.(49), and we can have 


nw = - 


273VT 3 4 5 

- -==r -|—r 2 + Cfc 

15^/777 5 


15 y/rj^ 


(r 3 + 2 r%) 


273 (arctan 1 ~ 2 ^ _ arcta n l + 2 ^ 




73 


73 


In 


(\/r + s/n) l (r + yfiV + r H ) 2 (r 2 + rr H + r 2 H ) 


( 68 ) 


(69) 


The integral constant Cu is hxed by the requirement that the final boundary stress tensor 
is in the Landau frame, which makes C k = -^T r J. So we finally have: 


k = [ F - i(r 3 + 2 r%)F ) 


(70) 
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In order to make a consistent check, one may put h, j and k into Eq.(58), it comes out 
that the three 1st order scalar perturbations that we have solved out satisfy Eq.(58) just 
right. So the scalar perturbations that we need to make Eq.(29) as the solution of Einstein 
equation again turns out to be 

ds(\)s — (—rdv 2 + rz Fh8ijdx l dx :j + 2r&Fjdvdr] dk/3k- (71) 

\r3 J 


D. Global form of the full metric containing first order perturbations 


Put all the stuff of zeroth and first order together, we get: 


r 4/3 


+ : 7§r) dv * + (dr 7/6 d v Pi - 2r n* dPdv 


+ 2r 1/6 (1 + FjdiPi ) dvdr + r 5/3 + -FSijdkfik + ) dx l dx j 

— 2 / r^x p, dp,f3 i dx l dr. 


(72) 


The above is just the full solution of the 1st order at the vicinity of aE = 0 in some special 
frame, whose covariant form can be constructed as: 


ds 2 = — r a ( f(r H (x),r) — 


Fk(r H (x),r ) 


d p u p I UpU u dx p dx u — 2rs(u ll a u + u u ap)dx p ’dx v 


+ rs F(r H {x),r)a^ u dx fJ, dx u + rs ^1 + ^F(r H (x), r)d p u p ) P fiu dx p, dx p 
— 2r^ (1 + Fj(r H (x),r)d p u p ) u^dx^dr, 


(73) 


where a pv = P^P^d^Uu) — \Ppvd p u p is the 4D covariant shear viscous tensor and a p = u v d v Up 
is the 4-acceleration relates with u„. 


IV. THE BOUNDARY STRESS TENSOR AND TRANSPORT PROPERTIES 

A. Derivation of boundary stress tensor 

The system of this model is a 5 dimensional Einstein gravity coupled with three scalar 
fields, its total action can be written as: 

S = Sb u ik + Sgh + S c .t., (74) 

where Sbuik is the bulk action (11), Sgh is the corresponding Gibbons-Hawking action: 

Sgh = —~2 [ d A xy/^-hK, 

J 


( 75 ) 
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with Hmn is the boundary metric tensor at a hyperplane with constant large r. K is the 
trace of the external curvature. The most crucial part in the total action is the bulk counter 
term S c _ t .. Since the bulk metric is not AdS, so the results of the counter term for AdS 
spacetime [72] cannot be directly used here, but fortunately there are also works on the 
renormalization of nonconformal branes [73]. Here we adapt the counter term used in [74], 
which they borrow from the much earlier work [28] on the renormalization of black D4 brane, 
in the Einstein frame it has the form of 

Sc.t. = \[ (76) 

«io J ^ 

Here H is the determinant of the boundary metric of Eq.(5) 

ds 2 = e~^' A h IJ , l ,dx^dx v + e 2A+SB dy 2 + e 2A ~ 2B dSt\. (77) 


Note x M = {a+,r}. After the dimensional reduction on the above metric, Eq.(76) becomes 


S c .t. = ~2 f d i x\T-h^- 

K 5 J ^ 


g 3^ 12 ^ 


This counter term contributes to the surface stress tensor as 


2 5S, 


C.t. 


y/^h dh> W K\ 


= —77 I —-e 12 


/IV 


(78) 


(79) 


Using Eq.(22) and remember that we’ve set L — 1, one has the surface stress tensor with 
contribution from the counter term as 


’ /IV 




T™ rf = A, ( K„ v - h,„,K - lr-^h 


/IV ,v /.LV J 


U /IV 


(80) 


In the standard technic for 3+1 decomposition of general relativity, h MN is defined as 


hMN — 9 mN ~ IT^mUn, 


(81) 


where Hm = NS7mt is the unit normal vector for a hyperplane at constant large r in the 5D 
bulk, of which the metric can be written as 


ds 2 = ( N 2 + N M N M )dr 2 + 2N M drdx M + h,MNdx M dx N . (82) 

N = {g AIN S7Nr)~^ is called the lapse function and N M is the shift vector. The index 
of hm and N M goes up and down with Hmn- The external curvature Kmn is related with 
h M N by 

Kmn = —-^dL n hMN = ( n p dpliMN + dMn p h PN + dNn p h PM ) , (83) 

in which C n is the Lie derivative along the unit normal n M . 
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B. Transport properties of QGP in D4 holographic QCD model 


The surface stress tensor that we obtain is 

1 


rj-isurf _ 

fils 


2 «! 


In 5 o 5 45 

r u Puu + ~r„u u u u - 2r 2 H cr^ - —r 2 H d p u p P pv 


2 2 

Comparing with the result in relativistic hydrodynamics: 


(84) 


- 2 ^ - oiyp,,,,. 


(85) 


where p, £, r/ and ( are the momentum density, the energy density, shear viscosity and bulk 
viscosity, respectively. We can get the respective hydrodynamical quantities for our system 
as 

1 1 3 1 5 3 1 | . 1 4 | 

P 2av§ 2 Hl " 2k\2 H ' 1 24 H: C 2ac| 15 H ' 

From Eq.(23) we can get the temperature for the 5D spacetime as 

1 

0~> 2 

T = ^X (87) 

As one can easily seen from the above two expressions, both the thermodynamic and the 
transport coefficients are only depend on temperature. This is due to the setup of this 
model. One can also get the entropy density as 

£ -\- V 1 - 

S = — = ^| 47Fr ^- ( 88 ) 

So the ratios of shear and bulk viscosity to entropy density are 



V = J_ C _ 1 
s 47t’ s 15n 


(89) 


Here we meet the renowned 1 /47r again and this suggests both the bulk of D4-D6 model and 
SS model belong to the class in [71], as [69] has pointed out. The bulk to shear ratio is 


c = ± 

77 15’ 


(90) 


which is also the same as in [69] and [70]. It is interesting to compare our result Eq.(90) 
with the results of Refs. [66, 67] of which for near-horizon non-extremal D4-brane is 1/10. 
This is understandable since the case we considered here is the compactified near-horizon, 
non-extremal D4-brane in which the relativistic fluid resides only on 1+3 dimensions out of 
the 1+4 dimensional D4-brane’s world-volume. The spacetime here comes from dimensional 
reduction on S 1 x S 4 . But in Refs.[66, 67], for D4-brane, the submanifold that is reduced is 
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the S ' 4 and the relevant hydrodynamics is 1+4 dimensional. Another consistency with [69] 
is the sound speed that can be obtained via thermodynamic quantities: 

o dp 1 

Cs= fe = 5' 

As a self-consistent check, we calculate the dispersion relations by using the constituent 
relation Eq.(84) as in [43]. If considering the temperature rn(x) and 3-velocity fa(x) has 
fluctuations as 



r H (x) = r H + Sr H e~ iuv+its , fa(x) = 5 fae~ iuv+its , (92) 

one can get the relations of the fluctuations by putting the above equations into the EOM 
of boundary fluid, i.e., the conservation equation for 

= 0. (93) 


Treating Sr h and 5fa as first order quantities, one can get the linear equation for the fluc¬ 
tuations 


5 

-ujSt h - r H kiSfa = 0, (94) 

y kiSr H + (r]( 2 k 2 - 3 ir H uj)Sfa + ^r^fakjSfa = 0 . (95) 

In order to make the above equations have non-trivial solution, the determinant of coefficients 
should be 0 , which gives 

c o =-shear mode 

3 + 

1 „ 4 

uj = ±—=\k\ — i -rTrP + 0(\k\ 3 ). sound mode (96) 

V5 15r^ /2 

Comparing with the following results in hydrodynamics 


fj 

u> = —i - k 2 , shear mode 

£ +p 


to = c s \k\ — i— — k 2 . 

1 ' 2 (e + p) 

sound mode 

(97) 

One can read the following relations 



V _ 1 c 2 _ 1 

C + tv 4 

(98) 

e + P 3+ 2 ’ ' 5’ 

2 {e + p) 15 r]i 2 ' 


Comparing with the results in Eq.( 86 ), we can find perfect consistency. 
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V. DISCUSSIONS AND OUTLOOKS 

We develop the BDE formalism of fluid/gravity correspondence in compactified black 
D4-brane background and investigate the transport properties of its gauge-side dual gluonic 
matter. Compactified D4 branes are background of the D4-D6 and the SS models, which are 
the two nonconformal top-down holographic QCD models. The SS model is a holographic 
model whose dual held theory lives on the world-volume of the flavor D8 branes, it is 
convenient to extract hadronic properties such as the meson and baryon spectrums from the 
SS model, since the DBI action of D8 branes describes some meson effective theory like the 
%-PT. However, in the SS model, people focus more on the flavor sector, and may ignore 
the bulk sector. Our current work focuses on the bulk sector, i.e., the compactified black D4 
brane background whose asymptotic region is not an AdS spacetime. Therefore, we choose 
the compactified black D4 brane background to describe nonconformal gluonic matter. 

The strategy is to use the dimensional reduction technique on the compact structure of 
SS model background in Einstein frame, and one can get a 5D effective Einstein gravity 
minimally coupled with 3 scalars with exponential potentials. Following the standard BDE 
formalism of fluid/gravity correspondence, we derive the constituent relation and read the 
thermodynamical and hydrodynamical quantities such as the energy and momentum density 
and the shear and bulk viscosities. It is found that the ratio of bulk to shear viscosity and 
sound speed from our results are consistent with the previous studies on the transport prop¬ 
erties of SS model [69, 70], which shows the validity of the BDE formalism of fluid/gravity 
correspondence in nonconformal background. The calculation of second order transport co¬ 
efficients are technically direct based on this work. What’s more, this work offers us an 
nonconformal prototype in fluid/gravity duality that is in parallel with Bharttacharyya et 
al.’s AdS^ construction, which provides us the opportunity to study nonconformal systems 
by the 5D metric (23) given in this work. 

As further applications, the most straightforward project is to calculate the second order 
transport coefficients of nonconformal gluonic matter by using the BDE formalism. The 
second order transport coefficients have been calculated in conformal systems with [46, 47] 
or without [43, 45] chemical potentials using the BDE formalism of fluid/gravity duality, 
and also been calculated via Green-Kubo formulae in both conformal [41] and nonconformal 
[75] systems. For the nonconformal Dp-brane backgrounds, the form of second order viscous 
tensor has been predicted in [67]. 

One may also study the effective fluid on a cavity with hnite r = R using the compactified 
black D4 brane solution (not under the near horizon limit) like [53, 54]. The most different 
point of the effective recipe from the present work is the dilaton, A and B should be boundary 
dependent since they both relate with // 4 , and thus rg 4 . This is like the case in [54] where 
scalars also relate with the harmonic functions, but unlike that in [53], in which scalar is 
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only r-dependent and the cut-off surface can be chosen as iso-dilatonic. Thus we should let 
those 3 scalars all have 1st order perturbations just like in [54], which may make us solve 6 
equations in all for the scalar part perturbation. 

Another interesting attempt in the future is to investigate nonconformal fluid with an 
axial chemical potential /x 5 via the SS model with smeared DO charge on the D4 brane world- 
volume [76], this model can extract the axial chemical potential and the axial charge diffusion 
constant besides the hydrodynamical quantities in the present work. What’s more, all the 
hydro dynamical quantities should be both temperature and the axial chemical potential 
dependent. We can also use the method developed in this work to investigate the newly 
found anomalous effects [77] (such as the Chiral Magnetic Effect, Chiral Separation Effect 
and so on) analytically. 
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Appendix A: Dimensional reduction from 10D to 5D 

We will do dimensional reduction on the following 10D diagonal background: 

ds 2 = e 2aiA g MN dx M dx N + e 2aaA {e 2plB dy 2 + e 2 ^ B -f ab d0 a d0 b ) . (Al) 

The non-trivial Christoffel symbol of this metric are (the ones with “tilde” are 10 dimensional 
components): 

f % P = T™ P + a^dpA + 5 AI d N A - g NP V M A)- 
fJJ = ~(a 2 V M A +/3 1 V M B)e^ 2ai+2a2)A+2 ^ lB - 
f V M y = ot 2 d M A + fadirfB; 

= -(a 2 V M A +p 2 V M B)e(~ 2ai+2a2)A+2 ^ B 7a6 ; 
f Mb = ( a 2$ mA + / 3 2 dMB)8 %; 

it = r a bc . (A2) 

From the above results, we also have: 

f mn = r mn + 5a:i<9jvfA, 

f mp + r y My + f Ma = + (5oq + 5 a 2 )dMA + (/31 + Af3 2 )dMB, (A3) 
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which can make our computation more convenient. The components of Ricci tensors are: 

Rmn =Rmn ~ (3«i + 5« 2 )V m V wj 4 — (/?! + 4/3 2 )VmVatR — ai5 , MivVpV p A 
+ (3 a 1 2 3 4 5 + 10aia 2 — 5al)dMAd^A — (3a 2 + 5aia 2 )gAfAr(ck4) 2 
+ (oti/3i — a 2 /?i + 4ai/3 2 — 4:a2/3 2 ){dMAdNB + OnAOmB) 

— {cxi/3i + Aai^2)gMNdpAd p B — ( f3 2 + 4 ^OmBOnB; (A4) 

Ryy = ~ \ a 2 V 2 A + f3{\7 2 B + (3aia 2 + 5a 2 )(<9A) 2 + (3ai/?i + 6a 2/ 5i + 4a 2 /3 2 )<9A ■ dB 
+ (/3 2 + 4 ^ 2 ){dB) 2 ] e (-2ai+2a a )A+2AB. (A5) 

= — [)a 2 V 2 A + /3 2 V'R + (3aia 2 + 50 : 2 )(<9A) 2 + (3o: 1 /3 2 T a 2 /?i + 9a 2 /3 2 )cL4 ■ dB 

+ (Plfa + 4/5|)(6>^) 2 ] e (-2a 1 +2a 2 )A+2 ft i ?7Q6 + 3^. (A6) 

Again, the components with “tilde” are 10 dimensional ones. Then we have the Ricci scalar: 

R m = [R~ (8ai + 10a 2 )V 2 A - (2/A + 8/A)V 2 R - (12a 2 + 30aia 2 + 30a 2 ) (dA) 2 
-6(a! + 2a 2 )(/A + 4 fc)dA ■ dB - (2/3 2 + 8/A/A + 20/3%) (dB) 2 ] e~ 2aiA 
+ i2e~ 2a2A ~ 2 ^ B . (A7) 
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